Shape oscillations of a charged diamagnetically-levitated droplet 
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We investigate the effect of electrical charge on the normal mode frequencies of electrically-charged diamag- 
netically levitated water droplets with radii 4.5-7.5 mm using diamagnetic levitation. This technique allows 
us to levitate almost spherical droplets and therefore to directly compare the measured vibrational frequencies 
of the first seven modes of the charged droplet with theoretical values calculated by Lord Raylcigh, for which 
we find good agreement. 



Understanding how electrically-charged liquid droplets 
vibrate, distort and shatter is important in high tech- 
nology areas such as ink jet printing and electrospray 
ionization for mass-spectromctry of biomolcculcs. The 
quiescent shapes and dynamics of the charged drop in- 
fluence the echos returned from rainfall radar, affect rain- 
drop breakup, and possibly lightning activity 1-3 . Based 
on some recent reports 4 ' 5 , we anticipate that the ef- 
fects of charge on the vibrations of 'liquid marbles', 
microliter powder-coated droplets with applications in 
microfluidics 6 , may also become an interesting topic for 
future study. The charged liquid droplet is the basis of 
Bohr and Wheeler's liquid drop model of the atomic nu- 
cleus and one might hope to gain new insights into the 
nucleus and its fission by studying the dynamics of this 
classical system. 

The starting point for understanding these areas is 
Lord Rayleigh's calculation of the frequency of the nor- 
mal modes of an electrically-charged droplet, made in 
1882 7 . However, it has been surprisingly difficult to ver- 
ify experimentally Rayleigh's result directly, owing to the 
requirement that the droplet be spherical and the mea- 
surement taken without contacting the liquid. 

Brazier-Smith et al. measured the oscillations of a sin- 
gle 1.05 mm radius water droplet suspended on a wire, as 
a function of voltage 1 , and Saunders and Wang verified 
this result for droplet radii up to 2.8 mm using air-flow 
to suspend the droplets 8 . They obtained good agreement 
with theory for the fundamental mode, but were unable 
to measure the oscillations of higher order modes. Trinh 
et al. obtained measurements of the first three modes 
using a hybrid acoustic and electrostatic levitation tech- 
nique, but were unable to compare their results directly 
with Rayleigh's theory, owing to the non-spherical equi- 
librium shape of the droplet levitated in this way 9 . The 
equilibrium shape of an acoustically-levitated or air-flow 
levitated droplet is non-spherical because the gravita- 
tional force on the whole droplet is balanced only at the 
liquid's surface. Hence, the droplet becomes distorted 
under its own weight. The frequency of the fundamental 
mode of small, 1-50 yU m diameter, nearly spherical wa- 
ter droplets, levitated electrostatically, has been shown 
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to agree with Rayleigh's theory 10-13 , but measurements 
of the frequencies of higher order modes have not been 
made using this technique. 

Here, we use diamagnetic levitation to study the nor- 
mal mode oscillations of an electrically-charged water 
droplet, with radii a between 4.5 and 7.5 mm, and surface 
electric charge of 10 -9 C <~ 10 10 e. Previous experiments 
have used diamagnetic levitation to study the dynam- 
ics of uncharged droplets 14-20 and formation of liquid 
films 21,22 . In contrast to acoustic, electrostatic or air 
flow methods, diamagnetic levitation balances the force 
of gravity at the molecular level throughout the body of 
the levitating fluid, effectively allowing one to study fluid 
dynamics in a nearly zero gravity environment. With this 
technique, one can levitate large mm to cm-scale droplets 
that arc nearly spherical at rest. The large size of the 
droplets allows us to make accurate measurements of the 
frequencies of the normal modes using a relatively sim- 
ple technique, and because the drops are spherical, we 
can compare the measurements directly with Rayleigh's 
theory. We observe the vibrations of the fundamental 
and higher order modes, up to the first seven modes and 
find good agreement with Lord Rayleigh's theory for all 
these modes. Whitaker et al. and Vicente et al. 23-25 
have measured the resonant frequencies of the first 14 
modes of diamagnetically levitated helium droplets, 1- 
3 mm in radius, with charge ~ 10 3 — 10 times smaller 
than in our experiments. However, in their experiments, 
the small charge was applied as part of a technique to 
excite vibrations by an electric field, and was too small 
to give a measurable charge-induced shift in frequency. 

We use a superconducting solenoid magnet with a ver- 
tical, room-temperature, 50 mm diameter bore to levi- 
tate the water droplets. A water droplet levitates ap- 
proximately 80 mm above the geometric center of our 
solenoid, where the diamagnetic force, proportional to 
BVB, is equal in magnitude to, and opposite in direc- 
tion to, the gravitational force on (i.e. the weight of) 
the droplet 26 ' 27 . The magnetic field is B w 12 T and the 
vertical field gradient is dB/dz w 120 Tm" 1 at the levita- 
tion point. We fine-tune the current in the solenoid until 
the equilibrium shape of the levitating droplet is as close 
to spherical as possible. This is achieved by reducing the 
quadrupolc component of the magnetogravitational po- 
tential 'trap' to zero; details of this procedure are given 
in Ref. 19. Droplets of the sizes we use here are spherical 
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FIG. 1. Cutaway drawing of a water droplet levitating within 
a spherical brass cavity inside a plastic sleeve within the 
50mm-diameter magnet bore. See text. 



action after the experiment. By measuring the difference 
in weight between the wet and dry paper, we obtained a 
second measurement of the droplet volume. Using this 
simple technique, we were able to determine that the 
mass loss through evaporation of the drop, during the 
measurement period (approximately 30-60 minutes) , was 
less than 2% for all the droplets we measured, except for 
the smallest (4.5 mm radius droplet), where the loss was 
3%. 

The procedure for charging the droplet is as follows. 
A thin wire, connected to earth, is brought into contact 
with the droplet through the hole in the top of the cavity. 
A voltage is applied between earth and the cavity enclo- 
sure, charging the droplet capacitively. Since the cavity 
is spherical, the charge on the droplet, as a function of 
voltage, is given by 



to a good approximation. We measure the ratio of the 
equatorial (horizontal) diameter to the polar (vertical) 
diameter to be 1.00 ± 0.02. 

Fig. 1 shows a schematic diagram of the water droplet 
levitating in the vertical magnet bore. The droplet levi- 
tates within a spherical metal cavity, which is electrically 
insulated from the magnet bore by a plastic sleeve. A 
voltage is applied between the cavity and earth by using 
a high voltage power supply. 

A 1 mm-diameter nozzle passes through the hole at 
the base of the cavity, and is directed at the center of the 
underside of the droplet, as shown in Fig. 1. It is con- 
nected to a rubber bulb outside the magnet by a tube. 
When the bulb is struck on a hard surface, the resulting 
pulse of air from the nozzle excites several shape oscilla- 
tion modes simultaneously, with amplitude less than 5% 
of the droplet radius. 

An optical fiber, inserted through the air nozzle, di- 
rects light from a HeNe laser at the droplet. The drop 
acts as a lens, focusing the light onto the aperture of a 
second fiber, through a hole at the top of the spherical 
cavity, as shown in the figure. This fiber transmits the 
light to a photodiode outside the magnet. The emf of the 
photodiode is recorded by a digital storage oscilloscope. 
Since the focal length of the drop depends on its shape, 
the intensity of the laser light falling on the photodiode 
oscillates as the drop vibrates. Other groups have em- 
ployed variations of this optical technique for measuring 
the shape oscillations (e.g. Refs. 11 and 24). 

The temperature of the water was brought to 16°C 
in a water bath before it was injected into the potential 
trap, to match the ambient temperature in the magnet 
bore. At this temperature, the density of the liquid is 
p — 999 kgm~ 3 , its surface tension is T — 73.3 mNm" 1 
and its kinematic viscosity is v = 1.11 x 10~ 6 m 2 s~ 128 . 
The liquid was injected into the trap using a glass pipette. 
The volume of liquid injected was determined to better 
than 1% uncertainty from the difference in the weight of 
the pipette before and after injection. The liquid was 
drawn out of the bore using a paper towel by capillary 



where R is the internal radius of the cavity, a is the ra- 
dius of the spherical droplet at rest, and V is the voltage 
between the wire and the cavity. 

The wire is slowly removed from the droplet and then 
the voltage on the cavity is returned to ground. The vol- 
ume of liquid withdrawn with the fine wire is negligible. 
Other methods have also been used to charge levitat- 
ing droplets (e.g. Ref. 29). The charge on the droplet 
leaks away with a time constant of several hours, allow- 
ing ample time to measure the normal mode frequen- 
cies. The maximum voltage applied between cavity and 
droplet was 3kV. Larger voltages caused the mechanical 
equilibrium between diamagnetic, electric and gravita- 
tional forces to become unstable, causing the droplet to 
escape the magnetogravitational potential 'trap' and hit 
the walls of the cavity. 

Fig. 2 top shows the oscillations in the photodiode 
signal developed by the light refracted through an a = 
7.5 mm droplet; they decay exponentially as exp(— \t) 
with a time constant 1/A = 0.1 — 10 s, dependent on 
a and the oscillation mode. The center panel of Fig. 2 
shows the power spectrum of these oscillations. Several 
peaks are evident in this spectrum, corresponding to the 
I = 2 to 5 Rayleigh modes. The peak near 2 Hz is due to 
a small oscillation of amplitude ~ 1 mm of the droplet 's 
center of mass of about the equilibrium levitation posi- 
tion. 

The lower panels show the power spectrum in three 
narrow frequency windows, centered on the I = 2, 3 and 
4 resonances. The resonances move to lower frequency 
with increasing charge Q. 

For small amplitude oscillations, the resonant frequen- 
cies are expected to be uj 2 — uiq + ljq, where 

_ Tl(l -1KT+ 2) 
W ° ^ (2) 

is the resonant frequency of the Zth mode of an uncharged 
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FIG. 2. Top: Oscillations in laser light intensity falling on 
the photodiode reveal the oscillations of a water droplet with 
radius 7.5mm. Center: The peaks in the power spectrum 
of the oscillations show clearly the 1 = 2 — 5 resonant fre- 
quencies of the droplet. The measured resonant frequencies 
are slightly higher than the frequencies predicted by Eqn. 2, 
indicated by dashed lines, owing to the effect of the magne- 
togravitational potential trap. Bottom: The peaks shift to 
lower frequency with increasing charge. Black (open circle) 
uncharged; red (filled circle) 1.4nC; green (square) 2.2nC; blue 
(triangle) 2.9nC; pink (cross) 3.6nC. 



droplet with surface tension T and density p, and 

Q 2 l(l- 1) 



167r 2 eopa 6 



(3) 



is the shift in frequency owing to the charge Q on the 
surface 7 . These equations assume that the viscosity is 
zero. We justify this assumption shortly. Since the 
Coulomb forces between the charges at the surface op- 
pose the cohesive force of surface tension, increasing 
the surface charge density reduces progressively the fre- 
quencies of the normal modes until the critical charge 
Q 2 R — 647r 2 eoa 3 T is reached, known as the Rayleigh 
limit. At this point the frequency of the fundamen- 
tal mode (I — 2) is reduced to zero, and the droplet 
becomes unstable to break up under the action of the 
Coulomb force 10-13 . Here, the largest voltage that we ap- 
ply charges the droplets to approximately Qr/3, which 
is approximately 1 x 10 10 e for a — 4 mm, and 3 x 10 10 e 
for a = 8 mm. 

The dotted vertical lines in Fig. 2 show the reso- 
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FIG. 3. Shift in frequency loq = oj(0)' 2 — oj(Q) 2 as function of 
charge Q for the I = 2 resonances. 



nant frequencies of an uncharged droplet, as given by 
Eqn. 2. The measured frequency of each resonance is 
slightly higher than we expect from Eqn. 2. This is 
due to the effect of the magnetogravitational potential 
trap, which acts as an additional cohesive force on the 
droplet. The magnetogravitational potential trap, within 
which the droplet levtitates, has the form U (r, 9) — 
co(r) + C2(r)P 2 (cos6') + C3(r)P 3 (cos 9), where r and 9 
are spherical polar coordinates centered on the levitation 
point {9 is the inclination) and P n are Legendre func- 
tions. In these experiments, the quadrupole component 
c-iPi has been reduced to zero by adjusting the current 
in the solenoid. There is no effect to first order in oscilla- 
tion amplitude and C3 on the frequencies. The spherical 
component of U adds a term c' (a)l/a to Eqn. 2, where 
c' (a) = dco/dr. This effect is discussed in Ref. 19 and 
20; see also Ref. 25. 

Fig. 3 shows the measured shift in frequency ufj — 
w(0) 2 — <^(Q) 2 as function of charge, for the I = 2 reso- 
nances, and Fig. 4, top, shows the shift in frequency 
for the I — 2 — 4 resonances. In the latter figure, we 
have plotted the dimensionless frequency versus the di- 
mensionless charge by scaling the frequency by the time 
scale t = (pa 3 /T) 1 / 2 and the charge by Q R . When the 
data are scaled this way, the resonances from droplets 
with different radii, ranging from a = 4.5 to a — 7.5 mm 
collapse onto the straight lines l>SqT 2 = 41(1 — 1)Q 2 /Q 2 R 
in agreement with Eqn. 3. 

When the frequency shift ujq is further scaled by 1(1 — 
1), as in the bottom panel of Fig. 4, the data for all the 
modes collapse onto the same line. This panel includes 
the measured peak positions of modes up to I = 8, in 
addition to the data shown in the top panel. 

The most significant contribution to the error in lj 2 -, 
comes from the uncertainty in the position of the peaks 
in the power spectrum. Although the uncertainty in the 
frequency of a particular resonance is approximately con- 
stant 5u> ~ 0.03a;, the uncertainty in the frequency dif- 



4 



mode / 



Frequency shift 
sealed by time scale r 2 = pa?/T 



O 2 
A 3 
□ 4 



0.04 







□ 


•tr' 

i □ 




] 


□ 

A 








0.06 


o.os 



Charge Q 2 /Q 2 , 



Frequency shift scaled by T 2 /l(l — 1) 




0.04 0.06 

Charge Q 2 /Q| 



FIG. 4. Top: Fan plot showing the shift in frequency lo c = 
uj(0) 2 — uj(Q) 2 with charge Q for the I — 2 — 4 resonances. The 
frequency has been scaled by the time scale r = (pa 3 /T) 1 ^ 2 
and the charge by Q 2 R — 64-7r 2 eoa 3 T. Bottom) when the fre- 
quency shift ojq is further scaled by 1(1 — 1), the data for all 
the modes collapse onto the same line. This panel includes 
the measured peak positions of modes I > 4. Only selected 
data points shown with error bars, for clarity. 



ference scales with to: 5uiq = 2(8ui/ui)(uj(Q) 2 +uj(Q) 2 ) ps 
A(Su) /uj)iJq. Hence the uncertainty in uiq is approxi- 
mately proportional to 1(1 — l)(l + 2), as indicated by 
the error bars on Fig. 4. 

The uncertainty in the volume of the droplet comes 
primarily from the evaporative loss, which is less than 2% 
for all the droplets we measured, except for the smallest 
(a — 4.5 mm), where the loss was 3%. We estimate the 
error in Q to be approximately 3% from the combined 
uncertainties in a and voltage. 

Although a strong magnetic field is expected to damp 
the motion of a conducting liquid through the field via 
the Lorentz force 30 , the conductivity of the distilled water 
used here is too small for any effect on the oscillation 
frequencies, or damping of the oscillations to be observed. 

The time constant A of the decay of the oscillations 
depends on the charge on the droplet. For example, for 
an a = 6 mm droplet, the I = 2 mode of the uncharged 
droplet has a time constant A = 0.16 ± 0.02 s _1 in good 
agreement with the theory A « vaT 2 (l — 1) (2Z + 1) for 
low viscosity liquids 31 . When we charge the droplet to 
(5^/3, the oscillations decay more quickly, with a time 
constant A = 0.26 ± 0.02 s _1 . We obtained the time con- 



stant for an individual mode by measuring the decay of 
the resonance peak in the power spectrum as a function 
of time. The faster decay of the oscillations of a charged 
droplet could be due to the action of the Lorentz force 
on the oscillating surface charges in the strong magnetic 
field. A similar effect has been observed in the oscilla- 
tions of charged helium droplets 24 . The Ohnesorge num- 
ber, which characterizes the ratio of viscous stresses to 
surface tension forces, is much less than unity for the sizes 
of water drops in these experiments. Hence, we expect 
the viscous damping to have an insignificant effect on the 
frequencies. For example, taking into account the viscos- 
ity of water, we expect the I = 2 mode of an a = 5 mm 
droplet to be just 5 x 10~ 4 % lower than that given by 
Rayleigh's formula 19 , Eqn. 2. Since the time constant of 
Lorentz- force damping is similar to that of viscous damp- 
ing, it is reasonable to neglect the action of the Lorentz 
force on the resonant frequencies as a small effect, too. 

In conclusion, we have demonstrated a method to 
study the dynamics of an electrically charged droplet us- 
ing diamagnetic levitation. Within the droplet, the force 
of gravity is compensated at the molecular level, such 
that the effective gravity acting at all points within the 
liquid is reduced to within a few percent of g. Because 
of this, diamagnetic levitation allows the levitation of 
nearly spherical droplets. Using this technique we have 
confirmed Lord Rayleigh's predicted effect of charge on 
the normal mode frequencies of a spherical droplet, for 
the first seven modes. 
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